Entanglement production in chaotic quantum dots subject to spin-orbit coupling 
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We study numerically the production of orbital and spin entangled states in chaotic quantum 
dots for non-interacting electrons. The introduction of spin-orbit coupling permit us to identify 
signatures of time-reversal symmetry correlations in the entanglement production previously un- 
noticed, resembling weak-(anti)localization quantum corrections to the conductance. We find the 
entanglement to be strongly dependent on spin-orbit coupling, showing universal features for broken 
time-reversal and spin-rotation symmetries. 
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I. INTRODUCTION 

The existence of entangled many-particle quantum 
states subject to nonclassical correlations is widely 
recognized as a fundamental resource for quantum 
information processing. Two quantum systems A and 
B are said to be entangled if they are not separable, 
namely, if the common (pure) state I^^ab) can not 

[ be written as the product of individual states I^^Ia) 
and \^b) (i.e., \^ab) ^ |4'a)|*b))- In such case the 

[ internal degrees of freedom of systems A and B are 
quantum mechanically correlated, since any measure- 
ment performed on system A would condition the results 
of a measurement on system B beyond any classical 
constraint. In solid-state physics, several aspects related 
to the production, control and detection of entangled 
electronic states have been addressed (see Ref. |l] for 

[ a recent review) and an extended literature already 

] ^..-^.l.^ 2.3.4.5.6.7.8.9.10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25 

Most of the existing proposals for the generation of 
electronic entanglement are based on the presence of 
some kind of interaction between the particles involved 
as, e.g., Coulomb interaction, (anti)ferromagnetism, 
superconducting pairing, etc. In contrast to the original 

] belief, it was recently recognized that interactions are 
indeed not necessary to produce entanglement. Non- 
interacting electrons, initially in a separable uncorrelated 
state, can evolve into an entangled state due to exchange 
correlations in a scattering process from an external 
potentiali2i2i^iSiSiL& This applies to electronic transport 
in multiterminal mesoscopic quantum conductors. In 
this approach, the efficiency of the entangler depends 
on the particular characteristics of the scatterer, which 

■ is described by its corresponding scattering matrix 
S. Among the several possibilities, systems of special 
interest are disorder-free chaotic quantum dots, also 
referred to as chaotic billards. These have the advantage 
of allowing for a statistical analysis that can reveal 
universal properties of (chaotic) electronic entanglers. 
This case was recently addressed by Beenakker et al^ 
by means of a random-matrix-theory (RMT) approach. 
They obtained a universal mean value for the degree of 
two-electron entanglement produced between spatially 



separated orbital channels, and remarkably found that 
it is not significatively affected by the breaking of 
time-reversal symmetry (TRS)- in contrast to other 
(single-particle) transport characteristics as the conduc- 
tance. Later, Samuelsson, Sukhorukov, and Biittikes^ 
reformulated the problem by formally including the spin, 
though in the absence of any spin-dependent interaction. 

In this paper we study the production of both orbital 
and spin two-particle entanglement for non-interacting 
electrons in multiterminal chaotic billards. We approach 
the problem numerically. By introducing the Rashba 
spin-orbit (SO) coupling^^ together with a magnetic flux 
breaking TRS, we are able to identify the signatures of 
weak-localization (WL) and weak-antilocalization (WA) 
quantum corrections in the entanglement production. 
Such TRS-correlation effects appear as a constraint that 
can limit/enhance the efficiency of the entangler. We 
find that the production of spin as well as orbital entan- 
glement is strongly affected by SO coupling on a scale 
corresponding to the pass from WL to WA. Additionally, 
we show that a finite residual entanglement survives af- 
ter the breaking of both TRS and spin-rotation symmetry 
showing some universal characteristics. 

The paper is organized as follows. In Sec. [H] we de- 
scribe our numerical model and shortly review some rel- 
evant single-particle properties as the WL and WA quan- 
tum corrections to the conductance. In Sec. IIIII we start 
discussing the interaction-free production of entangle- 
ment from a separable two-particle state (Scc. lIIlX|) . We 
introduce the concurrence as a measure of two-qubit en- 
tanglement in Sec. IIII CI In Sec. IIII Bl we discuss some 
features related to the accessible entanglement as a con- 
sequence of the local particle number conservation. The 
results are presented in Sec. lIVI followed by a short sum- 
mary of conclusions in Sec. 



II. MODEL AND SINGLE-PARTICLE 
PROPERTIES 

We consider a two-dimensional chaotic dot connected 
to electron reservoirs at the left and right, as shown in 
Fig.^S^i Two single-orbital-channel leads are attached at 
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FIG. 1: Chaotic quantum dot entangler used in the numerical 
model. Leads connected to electron reservoirs support one 
orbital channel plus two spin channels each. 



each side of the dot, similar to what proposed in Ref. % 
In addition, each orbital mode can support two spin chan- 
nels. The single-particle Hamiltonian for electrons with 
charge — e and effective mass m* reads 

where 11 = p + (e/c)A, or, is the Rashba SO couphng 
strength, cr is the vector of Pauli spin matrices, z is the 
axis perpendicular to the dot's plane, and 1^ is a hard- 
wall confining potential describing the dot's contour. The 
strength of the Rashba SO interaction can be given in 
terms of the ratio LescZ-^'SO) where Lose = -kA/w is the 
classical escape length in an open chaotic billard of area A 
with a total opening of width w,^^ and Lgo — ttIi^ / a^m* 
is the spin-precession length due to SO coupling. A uni- 
form magnetic field (generated by the vector potential 
A) introduces a flux (p (measured in units of the flux 
quantum 0o = hc/e). Applying a (small) bias voltage 
between reservoirs produces a coherent electron current 
through the dot from left to right. For the calculation of 
the corresponding scattering amplitudes we implement 
a recursive Green's function technique based on a spin- 
dependent tight-binding model arising from a real-space 
discretization.^^ Additionally, we also perform some in- 
dependent RMT simulations for comparison. 

The single-particle transport properties are character- 
ized by the Landauer-Biittiker linear conductance G. For 
illustration, in Fig.|21we present some numerical results 
for the sample-average conductance (G) of the chaotic 
dot of Fig. n as a function of Losc/Lso and (j). Coherent 
backscattering leads to a minimum in G (WL) at (/> = 
in the absence of SO coupling. For large SO coupling, G 
presents a maximum (WA) instead. The transition from 
WL to WA show s up around Losc/Lso ~ 10 (see also 
Refs. l30l3ll32l33j) . For a large 0, TRS is broken and G is 
independent of the SO coupling stregth, remaining close 
to its classical value (Gd — 2e^//i in our case). Such 
underlying single-particle physics is relevant for the un- 
derstanding of the new two-particle effects reported in 
this paper. The crossover from WL to WA also manifest, 
though differently, in the entanglement production as we 
see below. 
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FIG. 2: Sample-averaged conductance (G) for the chaotic 
quantum dot of Fig. vs. SO coupling in the presence of 
a magnetic flux (p. The curves correspond to 4'/4'Q = (solid 
line), (j)/(j)o = 1/2 (dashed line), and 4>/4>o ~ 5 (dotted line). 
The results illustrate WL (WA) quantum corrections to the 
classical conductance {2e^ /h in this case) for weak (large) SO 
coupling due to TRS {4'/4'o = 0). These quantum correction 
disappear as TRS breaks {4>/4>o ^ 0). 



III. ENTANGLEMENT PRODUCTION FROM 
SEPARABLE TWO-PARTICLE STATES 

A. Incoming and outgoing two-particle states 

We consider a separable two-particle state incoming 
from an electron reservoir on the left of the dot of Fig.^ 

|*in) = ar^ar^|0), (2) 

where a^'^ creates an incoming electron in lead i = 1, 2 
with spin Si =t, |, and |0) is the Fermi sea at zero tem- 
perature. Multiple scattering within the dot entangles 
the ongoing state. This is a coherent superposition of or- 
bital and spin channels determined by the single-particle 
^-matrix. It reads 

l^out) = EE^«i''^™2 K^b^^n. (3) 

n.a m,T 

where Sj^^ is the scattering amplitude from lead z = 1 , 2 
with spin Si to any lead j = 1,...,4 with spin s. The 
6^^ creates an outgoing electron in lead j with spin s, 
satisfying the vector equation ■ S = . The terms in 
(PI with n — m, (7 = T vanish for the sake of Fermionic 
statistics. 

B. Particle conservation 

The l^'out) of Eq. Q can be split into three terms with 
different local particle number at the left (til) and right 
(71r) of the dot such that til + = 2, in the form 

|*out) = ^ \nL,nn) 

= |2,0)-f |0,2) + |1,1). (4) 
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The accessible entanglemenfe^i in l^fout) is studied in Bell- 
like measurements by performing local operations that 
conserve the local particle numberii This is of fundamen- 
tal importance for electrons, preventing the local creation 
of states in a coherent superposition of different number 
of particles. This means that local operations do not mix 
the three terms of Eq. Q and the entanglement can be 
studied in each of them independently. 



C. The concurrence as an entanglement measure 

We evaluate the amount of entanglement between pairs 
of two-level (sub)systems or quhits, which in our case cor- 
respond to an electron leaving the quantum dot in one 
of two predetermined orbital/spin states. This requires a 
bipartition of the system by choosing some pairs of out- 
going channels of interest, tracing out the nonobserved 
degrees of freedom compatible with that choice. The 
two-qubit entanglement is quantified by the concurrence 
< C < 1, defined 



C{p) = max{0, Ai - A2 - A3 - A4}. 



(5) 



The AiS are the square roots of the eigenvalues (in de- 
creasing order) of the matrix pp, where p e 4 x 4 is a 
two-qubit density matrix and p = {ay ® <7y)p*{ay >S> (Jy), 
with ay the second Pauli matrix. Separable unentangled 
states have C ~ 0, while C = 1 correspond to maxi- 
mally entangled (Bell) states. States with < C < 1 
are non-separable, partly entangled states. A C ^ 
is a necessary and sufficient condition for affirming that 
the two involved qubits are non-classically correlated due 
to entanglement. For electrons in the quantum dot of 
Fig. where both orbital and spin degrees of freedom 
are involved, the two-particle density matrix p is gen- 
erally larger than 4x4, i.e., it does not correspond to 
a two-qubit system. Still, a study of the, e.g., orbital 
qubits can be performed by defining a reduced density 
matrix (RDM) T^sip) £4x4, where the trace is taken 
over the spin degree of freedom s; and vice versa (see 
Sec. lIVp . For a chaotic dot, the outgoing wave function 
(0) is a function of a random scattering matrix S, and so 
its corresponding density matrix. The entanglement con- 
tained in (PI depends on the particular S, which is sam- 
ple dependent. Hence, we characterize the production of 
entanglement in by calculating the sample-averaged (C) 
and its fluctuations var(C) = (C^) - (C)^. The con- 
currence can be related to zero-frequency current-noise 
measurements (i.e., without time- resolved detection) in 
mesoscopic conductors, as recently shown in Ref. 1 and 
U Alternatively, the presence of entanglement could be 
determined by other means as, e.g., beam-splitter cur- 
rent correlations giving a lower bound for entanglemenli^ 
or implementing some entanglement witness^ Here, we 
calculate C from its definition for several bipartitions 
of the outgoing state l|3Jl independently of any particular 
detection scheme. 




20 30 

Lesc/^SO 



FIG. 3: Left-right orbital entanglement: {C£r) vs SO cou- 
pling with TRS preserved (upper panel) and broken (lower 
panel). Solid (dashed) lines correrspond to (anti)parallel in- 
coming spins. Insets depict the fluctuations var(C), respec- 
tivelly. Full dots: RMT results from Ref. & 
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FIG. 4: Left-right spin entanglement: (Cl"^") vs. SO cou- 
pling with TRS preserved (upper panel) and broken (lower 
panel). Solid (dashed) lines correrspond to (anti)parallel in- 
coming spins. Fluctuations var(C) shown in insets. Full dots: 
RMT results. 



IV. RESULTS 



We investigate first the entanglement produced be- 
tween outgoing left and right channels. To this aim we 
project l^'out) onto the subspace containing one single 
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excitation at each side of the dot (i.e. |1, 1)). We obtain 

i^lr) ^ J2J2(^P^'^^J' - ^fr^^pT) ^p^&f io), (6) 

where the indices p = 1,2; a =t,i and q = 3,4;/3 =1" 
, J, stand for left and right outgoing channels, respec- 
tively. The density matrix of the state © reads plr = 
|*lr)(*lr|/(*lr|*lr) e 8 X 8, which by construc- 
tion accounts for both orbital and spin degrees of free- 
dom. The degree of orbital entanglement contained 
in jnj can be extracted from pLR by tracing out the 
spin degree od freedom, defining the RDM p'^^ — 
3('^'/'I/°lr|o!,/3) € 4 X 4 for two orbital qubits. The 
Fig.|3]shows results for the corresponding average concur- 
rence (ClI?) '^^ -^csc/-^so for parallel (si = S2 ~1; solid 
line) and antiparallel (si =1,S2 —i] dashed line) incom- 
ing spins subject to a magnetic flux 0. The insets depict 
the fluctuations, respectively. Outgoing spins of different 
species do not contribute to orbital entanglement.^ This 
is why only parallel incoming spins show a finite (C£5|'} « 
0.39 in the absence of SO couphng (iosc/-^so — 0) when 
leaving the quantum dot. This value, almost unaffected 
by the breaking of TRS (finite (j)), is in very good agree- 
ment with previous RMT results^ depicted by the full 
dots in Fig. 121 As SO coupling increases, spins flip dur- 
ing transport and outgoing spin channels of different sign 
open up both at the left and right of the dot. This hinders 
the production of orbital entanglement from originally 
parallel spins, leading to a reduction of the concurrence. 
In contrast, for incoming antiparallel spins SO scattering 
contributes to the formation of orbitally entangled states 
between spins of the same outgoing species. The scale on 
which the concurrence varies signiflcativelly as a function 
of Lesc/ Lso is similar to that determining the transition 
fron WL to WA in the conductance shown in Fig.|21 For 
large SO coupling, the degree of entanglement saturates 
as the orientation of outgoing spins randomize. However, 
a finite difference AC w 0.05 survives between different 
incoming spin configurations for (j) — (Fig. |31 upper 
panel). This is a consequence of the TRS correlations 
preserved by the SO interaction: As soon as a finite 
breaking TRS is applied the concurrence tends rapidly 
to a common asymptotic value independent of the ini- 
tial condition (Fig. O lower panel). This indicates that 
breaking time-reversal and spin-rotation symmetries give 
rise to a residual orbital entanglement with universal av- 
erage concurrence (Cj^^} w 0.075 for chaotic dots. We 
point out that orbital entanglement is very sensitive to 
the spin dynamics even for broken TRS (Fig. |31 lower 
panel), where WL and WA quantum corrections to the 
cond uctance are absent (dotted line in Fig. El see also 
Refs.lSO lSllslsl . Regarding the fluctuations (insets in 
Fig.jSJ, they show a functional dependence similar to that 
of the concurence. For large icsc/iso, \/var(C) « (C). 
These features repeat in our results of Figs. ^ and 
We further note that some difficulties may appear for 
detecting orbital entanglement produced from incoming 
antiparallel spins by violating Bell inequalities for shot 
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FIG. 5: Transmitted spin entanglement: {C'^'"'} vs. SO cou- 
pling with TRS preserved (upper panel) and broken (lower 
panel). Solid (dashed) lines correrspond to (anti)parallel in- 
coming spins. Insets show the fluctuations var(C). Full dots: 
RMT results. 

noise as described in Ref. 2. This is because both incom- 
ing as well as outgoing channels would mix at the left side 
of the dot (unless they can be spatially separated) . How- 
ever, this does not exclude the possibility implementing 
alternative approaches for detection as, e.g., the deter- 
mination of lower bounds for entanglement from beam- 
splitter current correlations^^ or the use of some entan- 
glement witness 1^ 

Information regarding the degree of spin entanglement 
between left and right channels contained in ^ can 
be evaluated from plr by tracing out the orbital de- 
gree of freedom instead, constructing the RDM Plr" ~ 
'zIplrIp, q) e 4x4 for two spin qubits. Results for 
the corresponding average concurrence (C^^") are pre- 
sented in Fig. 21 In contrast to the previous case of orbital 
entanglement, antiparallel incoming spins (dashed lines) 
lead now to a finite (C^^) already at Lcsc/ Lso = due 
to exchange correlations, while parallel incoming spins 
(solid lines) do not. The result is in agreement with our 
independent RMT simulations (full dots) performed by 
following the approach of Ref. [j. The presence of multi- 
ple orbital channels give rise to an outgoing mixed state 

(TrpLfi"^ < 1) with (Cl^") < 1. This differs from the 
case in which one single-orbital-channel lead is attached 
at each side of the dot: There, antiparallel incoming spins 
escape at the left and right in a pure singlet state with 
(^lr") ~ 1 independently of the scattering amplitudes 
(straightfoward from Eq. ©; see also Ref. 5). We also 

note that (Cl^") ^ (C'lr) Lcsc/ Lso = (compare 
Figs. 121 and This is probably related to the fact 

that, in contrast to spin entanglement, orbitally entan- 
gled electrons leave the dot at left and right in a pure 



5 



state (TrpLR = !)■ For large SO coupling, Fig. ^ shows 
features similar to those for orbital entanglement: A fi- 
nite AC survives between different incoming states at 
= due to TRS correlations (Fig. 0] upper panel). 
The difference disappears as TRS is broken by a finite 
(j) (Fig. 21 lower panel). More interestingly, the asypm- 
totic value for (C^^") is very similar to that for {Cf^"^) 
in Fig. 121 (lower panel). This indicates the existence of a 
universal value for the concurrence of residual left-right 
entanglement independently of the initial condition and 
particular degree of freedom. The claim is supported 
by RMT, which shows that when both time-reversal and 
spin-rotation symmetries are broken the ^-matrix is uni- 
formly distributed in the unitary group and no effective 
difference exists between spin and orbital channels. '^^ 

We consider now the entanglement production for 
transmitted spins, i.e., the spin entanglement between 
channels at the right side of the dot. This entangle- 
ment is contained in the |0, 2) component of Eq. We 
note that such component splits at the same time into 
other three contributions, each of them with different lo- 
cal particle number at leads 3 (n^) and 4 (n^) such that 
+ n4 = 2. Following Sec. 1111 lj| we see that also here 
the entanglement can be studied independently in each 
term due to local particle number conservation. The only 
component of interest is that one with one single excita- 
tion on each lead, obtained by projecting the outgoing 
state ^ onto the corresponding subspace. It reads 



|vI/r)^^(537^5: 

a,/3 



I3S2 

42 



^ff^^sT) fcg-t&f |0), (7) 



where a and P label ongoing spins in the leads 3 and 
4, respectively. Note that the |^r) can be only spin- 
entangled since there are just two orbital channels on the 
right side of the dot. After defining the density matrix 
^spm _ |^r)(5'r|/(4'r|4'r) for the two spin qubits, we 
plot in Fig. (upper panel) the corresponding {C^™) 
vs Lpsc/Lso for ip — 0. As in the previous case of left- 
right spin entanglement, only antiparallel incoming spins 
(dashed line) lead to a finite (C^™) at Losc/^so = 0. 
However, the degree of entanglement is much larger (its 
value agrees with our RMT calculations depicted by the 
full dot). That holds true for large SO coupling, where 
the concurrence arrives at a relatively large common 
asymptotic value (C^™) « 0.4 independently of the ini- 
tial condition. This contrast with our findings for left- 
right entanglement, where TRS correlations are relevant. 
The application of a finite (f), Fig. O (lower panel), does 
not affect the zero-flux characteristics significatively. 

Usefull information can be still extracted from the 
wave functions © and ((TJ by plotting their modulus 
square (Fig. 1^. The larger contribution is given by 
^'lrP- The I^'rP is much smaller instead. It means 
that highly spin-entangled electrons transmitted to the 
right are actually produced with a lower probability. This 
can be understood by using classical probabilities in the 
limit of broken time-reversal and spin-rotation symme- 
tries: The probability for the two incoming electrons to 
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FIG. 6: Sample-averaged probabilities for the states (|SJ and 
Q vs. SO coupling. Left (right) panel shows results for 
preserved (broken) TRS. Solid (dashed) lines correrspond to 
(anti)parallel incoming spins. 



end up in any two different channels is 1 /28. There are 16 
combinations giving one electron on the left and one on 
the right, and only 4 with one electron in each lead 3 and 
4. The quantities 16/28 « 0.571 and 4/28 « 0.143 are 
in agreement with the asymptotic values of Fig. El (lower 
panel) . Quantum corrections to these values show up for 
small SO coupling, specially when </> = 0. We further note 
in Fig. El that a probability difference appears between 
parallel (solid line) and antiparallel (dashed line) incom- 
ing spins in the left-right component due to Fermionic 
statistics. Breaking TRS increases the contribution keep- 
ing the relative difference unaffected. 



CONCLUSION 



In summary, we studied the role of spin dynamics and 
TRS correlations in the production of entanglement in 
mesoscopic conductors and the connection with WL and 
WA quantum corrections. By including SO coupling, 
among other things we found that the TRS effects can 
be more important than originally thought^ These man- 
ifest both in the degree of entanglement as well as in 
the production rate. The effects of SO coupling appear 
on a scale corresponding to the pass from the regime of 
WL to WA in the quantum conductance. We also deter- 
mined some universal characteristics of chaotic entanglers 
as the residual amount of entanglement produced after 
time-reversal and spin-rotation symmetry breaking. 
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